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Abstract 

The transition to intermittent mean-field dynamos is studied using numerical sim- 
ulations of isotropic magnetohydrodynamic turbulence driven by a helical flow. The 
low-Prandtl number regime is investigated by keeping the kinematic viscosity fixed 
while the magnetic diffusivity is varied. Just below the critical parameter value for 
the onset of dynamo action, a transient mean-field with low magnetic energy is ob- 
served. After the transition to a sustained dynamo, the system is shown to evolve 
through different types of intermittency until a large-scale coherent field with small- 
scale turbulent fluctuations is formed. Prior to this coherent field stage, a new type 
of intermittency is detected, where the magnetic field randomly alternates between 
phases of coherent and incoherent large-scale spatial structures. The relevance of these 
findings to the understanding of the physics of mean-field dynamo and the physical 
mechanisms behind intermittent behavior observed in stellar magnetic field variability 
are discussed. 

keywords: magnetic fields - MHD - turbulence. 

1 Introduction 

The observation of strong magnetic fields in astrophysical bodies (planets, stars, and 
galaxies) suggests the existence of a dynamo process, whereby a weak (seed) mag- 
netic field is amplified due to the conversion of kinetic energy into magnetic energy. 
Dynamos can be classified as large-scale (or mean-field) dynamos, or small-scale (or 
fluctuation) dynamos, according to whether the magnetic fields grow on spatial scales 
larger or smaller than the energy carrying scale of the fluid motion. A typical mani- 
festation of a large-scale dynamo is the solar cycle, where the distribution of sunspots 
in space and time display large-scale spatial coherence and long-term temporal corre- 
lation, as seen in the butterfly diagram 51, 52} [57] . Although the maxima and minima 
of solar activity form a recurrent 11-year cycle, long periods of very low solar activity, 
grand minima such as the Maunder minimum [3], have led several authors to look for 
a description of the solar cycle as an intermittent event due to the chaotic nature of 
the dynamo, driving the system to random alternations between phases of "regular 
magnetic activity" and grand minima. Many works employ low-dimensional models 
based on ordinary differential equations to investigate these modulations [e.g., Covas 



et al. [TB] and Wilmot-Smith et al. [BUj]. Others have focused on high-dimensional 
mean-field models based on partial differential equations. Covas & Tavakol 17 re- 
ported the existence of crisis-induced and Pomeau-Manneville type-I intermittencies 
[3111121 fHfl 114] in an axisymmetric mean-field dynamo model, suggesting the multiple- 
intermittency hypothesis, since more than one type of intermittency may be responsi- 
ble for the minima observed in solar cycle data. Ossendrijver [57] showed that a 2-D 
mean-field dynamo model that features an Q-effect based on the buoyancy instability 
of magnetic flux tubes could produce grand minima. The grand minima observed by 
this model were later described by Ossendrijver & Covas [3S] as a manifestation of 
crisis-induced intermittency [211 1421 1431 1461 147] . Intermittency in a 2-D mean-field 
model was also reported by Moss & Brooke [33] , where the reaction of the Lorenz 
force on the rotation is used as the nonlinear effect that limits the magnetic field 
at finite amplitude. Charbonneau [10] showed that an axisymmetric 2-D solar cycle 
model based on the Babcock-Leighton mechanism of poloidal field regeneration can 
exhibit intermittency in the presence of low-amplitude noise. Recently, Brandenburg 
& Spiegel [8] observed on-off intermittency in a mean-field dynamo model after im- 
posing stochastic fluctuations in the a-effect or using a fluctuating electromotive force. 
A review on the use of intermittent chaotic models to capture the main qualitative 
aspects of the temporal and spatial variability of the solar-cycle is provided by Spiegel 
153]. 

This paper employs 3-D magnetohydrodynamics (MHD) simulations to investigate 
the onset of intermittent mean-field dynamo as a function of the magnetic Prandtl 
number Pr, defined as the ratio between the kinematic viscosity v and the magnetic 
diffusivity r\. Although the model geometry is highly idealised, the full induction and 
momentum equations are solved, so that there is no need to appeal to an averaging 
process. Early numerical works in similar models focused on Pr > 1 regimes due to 
difficulties in exciting a nonlinear dynamo for low values of Pr [see, e.g., Nordlund et 
al. [36]]. Nevertheless, this regime is crucial for understanding astrophysical plasmas, 
where Pr is usually much less than one [Pr ~ 10 _7 -10~ 4 in the solar convective zone, 
depending on the depth |49|] . Furthermore, the role of compressibility of the veloc- 
ity field in astrophysical plasmas should not be ignored, even in low-Mach-number 
regimes [481 122] - Thus, we adopt a compressible MHD code and explore low-Pr 
regimes. Two conditions have been shown to be important for a fluid flow to act as 
a large-scale dynamo: chaotic stream-lines and kinetic helicity [111 [5] . Depending on 
the characteristics of the velocity field and other parameters, such as r), the resulting 
magnetic field lines may display regular or irregular motion, or a decay to a purely 
hydrodynamic state. In order to obtain a chaotic and helical velocity field, the flow 
is driven by an ABC (Arnold-Beltrami-Childress) forcing, which is a superposition of 
three helical waves [HE] with a characteristic wave number kf, which sets the energy 
carrying scale of the flow. The emergence of a mean field is studied by forcing the flow 
at scale kf = 5 and observing the energy transfer towards larger scales. 

This paper is divided as follows. Section [2] contains a description of the dynamo 
model. The main results are described in section [3] The kinetic viscosity is set at 
a value such that the velocity field is weakly turbulent. Then, a seed magnetic field 
is applied and the magnetic diffusivity r\ is progressively reduced until the onset of 
dynamo is observed. Three types of dynamo regimes are observed as a function of r\. 
In the first regime, right after the onset of dynamo, there is an intermittent switching 
between bursty phases of high-amplitude magnetic activity and quiescent phases of 
low magnetic energy, similar to the on-off spatiotemporal intermittency reported in 
other works [551 1561 1431 147] . The second regime reveals an intermittent switching 
between coherent and incoherent large-scale structures and, up to our knowledge, 
has not been reported in previous dynamo studies. This intermittent process persists 
up to a certain threshold value of the magnetic diffusivity, where the third regime is 
observed, with the system self-organizing into a spatially coherent mean-field that 



2 



exhibits complex temporal dynamics. An analysis of the spatiotemporal complexity 
of the patterns observed in each regime is also described. The conclusions are given 
in section [4] 



2 THE MODEL 

We consider a compressible isothermal gas (the ratio of specific heats 7 = c p /c v — 1, 
where c v is the heat capacity at constant pressure and c v is the heat capacity at con- 
stant volume) with constant sound speed c s , constant dynamical viscosity fi, constant 
magnetic diffusivity 77, and constant magnetic permeability po- The continuity equa- 
tion is solved in terms of the logarithm of the density In p, since this quantity varies 
spatially much less than density. 

^^+uVlnp + V-u = 0, (1) 
at 

where u is the fluid velocity. The momentum equation is given by 

1 u-Vu = ^ + + — I V u + - VV • u ) + f , (2) 



dt p p p V 3 

where J = V x B/po is the current density, p is the pressure and f is an external 
forcing. The pressure gradient is obtained from the entropy equation for an ideal gas 
s = c v ln(pp~ J ). Thus, Vp/p = c 2 lnp, where c 2 = 7p/p is assumed to be constant. 
The induction equation is written in terms of the magnetic vector potential A, so that 
V ■ B = 0, since B = V x A 

8A 

— = u x B - rjfj, J. (3) 
dt 

We adopt nondimensional units, such that c s — po — fio = 1, where po = (p) is 
the spatial average of p. Equations {l)-^ are solved with the PENCIL CODeJ^J in 
a box with sides L = 2n and periodic boundary conditions. The PENCIL CODE 
is a compressible MHD code that employs high-order finite differences and has been 
extensively used in astrophysical simulations [see Brandenburg & Subramanian [7] and 
references therein]. The initial condition is lnp = u = 0, and A is a set of normally 
distributed, uncorrelated random numbers with zero mean and standard deviation 
equal to 1CP 3 . For the forcing function f we adopt the form of ABC flow used in Sur 
et al. [54] 

, / sin kfz + cos kfy \ 
f(x) = — L I sin kfX + cos kfZ 1, (4) 
^ \ sin kfy + cos kj-x / 

where Af is the amplitude and kf the wavenumber of the forcing function, which is 
isotropic with respect to the three coordinate directions. The ABC forcing is an inter- 
esting choice for dynamo studies, since it is a Beltrami flow, Vu oc u, and, therefore, 
is maximally helical, i.e., (-fffc) 2 = (u ■ V x u} 2 = (|u| 2 ^ ^|V x u| 2 ^, where the angle 
brackets denote spatial integration over the periodic domain. 

In all the following sections we use Af = 0.1, which yields Mach numbers below 
0.5. A numerical resolution of 64 3 mesh points is chosen. As in Brandenburg [H], we 
set kf — 5 in order to be able to see the emergence of a large scale magnetic field, with 
spatial scales larger than the energy injection scale. Spatial averages are denoted by 
(•} and time averages by {•).. Unless otherwise stated, references to kinetic (Re) and 
magnetic (Rm) Reynolds numbers are based on the forcing scale 
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u x u y u z 

Figure 1: (Colour online) Contour plots of the velocity components for hydro- 
dynamic simulations at (a) v — 0.02 and (b) v — 0.005. 



Re=^L, R m =^, (5) 
v r/ 

where v = n/po is the average kinematic viscosity, A/ = 2n/kf is the forcing spatial 
scale, and U = (u 2 ^ 1 ^ 2 is the mean velocity at a time when the magnetic field is 
saturated. 

3 RESULTS 

3.1 Hydrodynamic simulations 

We start with hydrodynamic simulations in the absence of magnetic fields. We set 
kf — 5, A = 0.1, and vary the average kinematic viscosity v. Figure [TJa) displays the 
velocity components of an asymptotically stable solution of Eqs. (JT|-([2| at v = 0.02, 
which corresponds to Re ~ 12.38. The initial velocity field u = evolves with time 
until it takes the shape of the ABC flow. After converging to this attractor, the 
amplitudes of the velocity components do not change with time, so this is a steady 
state of the system. As the kinematic viscosity is reduced, the ABC flow becomes 
unstable and a sequence of symmetry breaking bifurcations takes place. At v = 0.005 
(iie ~ 100 based on the forcing scale, or Re ~ 500 based on the box scale), the 
system is in a weakly turbulent regime, as illustrated in Fig. [ljb) . The flow is weakly 
compressible, as can be seen in Fig. [2] which shows that the density fluctuations are 
within 5% of the mean value po — 1. 

3.2 Onset of dynamo action 

We now ki/ = 0.005 and change the magnetic diffusivity rj as we look for the onset of 
nonlinear dynamo action. Since we are interested in the low-Prandtl number limit, rj is 
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ln(p) 

Figure 2: (Colour online) Contour plot of ln(p) for the hydrodynamic simulation 
at v = 0.005. 

made larger than v, as Pr is varied from 0.0625 to 0.5. Figure [3] shows the bifurcation 
diagrams for the time-averaged magnetic ({E m ) t , red triangles) and kinetic ({Ek) t , 
black circles) energies in linear [Fig. |3|a)] and linear-log [Fig. [3jb)] scales as a function 
of r\ (lower axes) or Rm (upper axes). For each value of rj, the initial variations of 
B 2 ) /2 and (u 2 ) /2 are discarded, and time-averages are computed for the saturated 
energy values. For large values of rj, the seed magnetic field decays rapidly and there 
is no dynamo. At the onset of dynamo action at rj ~ 0.053 (Rm ~ 9.5 in the forcing 
scale, or Rm ~ 47.5 in the box scale), the magnetic energy starts to grow at the 
expense of kinetic energy. Two different dynamo behaviours can be identified in Fig. 
[3] In the first range, between rj ~ 0.053 and rj ~ 0.04 there is a sharp increase in 
the saturated magnetic energy. Between rj ~ 0.04 and rj ~ 0.01 the magnetic energy 
increases more slowly as an exponential function of rj until it is comparable to the 
kinetic energy at rj — 0.01. 

We turn to the dynamics near the transition to nonlinear dynamo. Figure [4] plots 
the time series of magnetic energy for different initial conditions at rj — 0.055, just 
before the onset of dynamo action. For a small seed magnetic field there is an initial 
exponential growth of the magnetic energy [Fig. |4ja)] , as expected from kinematic 
dynamo theory. Then the field saturates and starts to decay. Even during the growth 
phase the magnetic field has a much lower magnitude than the velocity field. Con- 
sequently, the impact of the Lorentz force on the velocity field is negligible and the 
same kind of behaviour was observed in kinematic simulations, where the Lorentz force 
term J x B/p was removed from Eq. |2|. Figures |4|b)-[4jd) show other instances of a 
"transient dynamo" for different types of initial conditions. Thus, the dynamo state 
is not an attractor of the system and we believe this transient and apparently chaotic 
dynamics displayed by the magnetic field is a signature of nonattracting chaotic sets 
known as chaotic saddles [26, 44, 45 14 , which have attracted wide attention recently 
due to their role in transition to hydrodynamic turbulence [391 1241 1591 125] . 

3.3 Intermittent dynamo 

At rj = 0.053, there is a transition to sustained dynamo action. The time series of 
magnetic energy exhibit random switching between phases of bursty and quiescent 
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Figure 3: (Colour online) Variation of the time-averaged kinetic (black circles) 
and magnetic (red triangles) energies as a function of r\ in (a) linear and (b) 
linear-log scales. The upper scales display the magnetic Reynolds number based 
on the forcing scale. 
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Figure 4: Time series of magnetic energy at rj 
dynamics for different initial conditions. 



0.055 exhibiting transient 



activity, as seen in Fig. |5|a), including long periods of minima with almost zero 
magnetic activity such as the one between t ~ 17500 and t ~ 20000. This is the same 
type of behaviour reported by Sweet et al. [551 156] for low Reynolds number (Re = 
6.3) dynamo simulations with Prandtl number close to one and ABC-forcing at scale 
kf — 1. A blowout bifurcation, whereby the hydrodynamic state looses transversal 
stability, was characterized as responsible for the intermittency, which in this case 
is called on-off intermittency, since the solutions go arbitrarily close to a manifold 
defined by the purely hydrodynamic state B = ("on" phases) and suddenly depart 
from the manifold during the strong magnetic bursts ("off" phases). 

For n = 0.05, the dynamo is strongly intermittent, but no grand minima are 
found in the time series of (B 2 ) /2, shown in Fig. ^b). There are strong bursts 
interspersed by lower peaks, but the energy is seldom close to zero. A look at the spatial 
structures of the magnetic field is helpful to distinguish the two types of intermittency 
exemplified by Figs. |5^a) and[5^b). Figure [6] shows the contour plots of B x for the 
on-off intermittency (rj = 0.053) at five different times. The snapshots in the left 
column (predominantly red) are taken at quiescent phases, and the ones in the right 
column, at bursty phases. The magnetic field displays complex spatial structures in 
both phases, although they are more clearly seen in the bursty phases. The dynamics 
at rj — 0.05 is very different. The contour plots of B x , shown in Fig. [7] reveal that 
a large scale sinusoidal modulation of the magnetic field is evolving, but there is an 
intermittent switching in the preferred direction selected by the magnetic field. In the 
first snapshot (t — 5000) the preferred direction for B x is z; in the second snapshot 
(t — 6000) the field is restructuring and there is no preferred direction; in the third 
snapshot (t = 9000) the preferred direction is y and there is no preferred direction at 
t = 10000. This coherent -incoherent intermittency has been observed for all tested 
values of rj in the interval from rj = 0.05 to r\ = 0.02. 

For rj = 0.01 the magnetic field has settled to a large scale and approximately 
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Figure 5: Intermittent time series of the magnetic energy at (a) r\ — 0.053 and 
(b) rj = 0.05. 



sinusoidal mean field, with complex small scale spatial structures and irregular oscil- 
lations, as shown in Figs. [8]and[9] In Fig. j^a) the magnetic energy (-B 2 ^ jl is strong 
enough to cause a reduction of the mean kinetic energy from about (it 2 ) /2 = 0.4 in 
the previous examples, to (w 2 ) /2 = 0.3. There is also a definite correlation between 
the time series of the velocity and magnetic field components. 

3.4 Inverse cascade 

The appearance of a mean-field dynamo is due to the transfer of magnetic energy 
from the energy injection scale to larger scales. This transfer is expected in helical 
flows and in mean-field dynamo theory is attributed to the a-effect [331 [29] , whereby 
interactions between low-scale fluctuations of B and u produce a large-scale B [other 
mechanisms can be responsible for the rise of a mean-field in the absence of a mean 
a-effect, such as the interaction of a fluctuating a-effect and large-scale shear flows 
[411 127] . The kinetic and magnetic energy power spectra for the coherent mean-field 
state at 77 = 0.01 are shown in Fig. [To] for t — 5000. The one-dimensional spectra 
are obtained from the three-dimensional spectra by computing the integrals of the 
spectral energy along spherical shells with rays defined by the modulus of k. For the 
magnetic spectrum (solid red line), although the forcing scale is kj = 5, there is a 
backward transfer of magnetic energy, making k — 1 the predominant scale, resulting 
in a B field with a typical scale of the size of the box. As for the kinetic energy 
(dashed black line), k = 5 is still the predominant scale. For this value of Rm the 
spectra are particularly similar in the inertial range (not properly defined yet, since 
they are not sufficiently extended). The k~ 5 line is a guide to the eye, and is the 
power-law observed for the magnetic spectrum of dynamos with helical flows at low 
Rm (following a k~ z range) .35, 132]. 

Figure [TTJplots the peak height for the Fourier modes k — 1 (dashed line) and k = 5 
(solid line) of the time-averaged magnetic energy spectra as a function of 77. It shows 
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Figure 6: (Colour online) Contour plots of B x for 77 = 0.053 at four different 
values of t. The magnetic energy increases and decreases randomly in time, as 
in Fig. |5ja) , resulting in on-off intermittency. 



9 



r| = 0.05 

t = 5000 t = 6000 




t = 9000 



5- 



0.05 



-0.05 




no. 05 




-0.05 




0.05 



-0.05 



t= 10000 




0.05 



-0.05 



Figure 7: (Colour online) Contour plots of B x for rj = 0.05 at four different 
values of t. There is an intermittent switching between large scale approximately 
sinusoidal states (at t = 5000 and t = 9000, for instance) and incoherent states. 
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Figure 8: (Colour online) (a) Time series of the magnetic (red) and kinetic 
(black) energies at ry = 0.01; (b)-(d) time series of the magnetic (red) and 
velocity (black) field components at a given point in the 64 3 numerical mesh. 



that the energy difference between these scales increases as the magnetic diffusivity 
decreases, giving rise to a progressively more coherent mean-field, as discussed in the 
next section. The two distinct ranges previously identified in Fig. [3] can be seen in 
Fig. [Tl]as well. The first one with a steepest energy increase between r] = 0.053 and 
r\ = 0.04, and the second one between r\ — 0.04 and r\ = 0.01, where energy seems to 



increase exponentially, as seen in the linear-log plot of Fig. 11 b) 



3.5 Spatiotemporal complexity 

Two operators computed in the Fourier space can be used to quantify the spatiotem- 
poral complexity of a system. The amount of spatial disorder can be quantified by 
means of the spectral entropy [401 1611 147| 

M 

S(t) = -£>*(*) MPk(*))> (6) 

k=l 

where M is the number of modes and Pk (t) is the relative weight of a Fourier mode 
at an instant t 

Mt) = v^nfw- (7) 

Since u and B are real variables, 

|6_ h (t)| = |6 fc (t)|, (8) 

and only Fourier modes with k > need to be considered. Mode bo(t) is decoupled 
from the other modes and is null for all t if &o(0) = 0. The spectral entropy is maximum 
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Figure 9: (Colour online) Contour plots of B x for rj = 0.01 at four different 
values of t. A robust large-scale coherent structure can be observed throughout 
the simulation. 



12 




Figure 10: (Colour online) Magnetic (solid red line) and kinetic (dashed black 
line) energy spectra at t = 5000 for 77 = 0.01. The fc~ 5 line is a guide to the eye. 



for a random system with uniform distribution, i.e., for all k, Pk(t) = 1/M, in which 
case S(t) = InM For the present spatial resolution, M = 32 and the maximum 
entropy is Sm ax ~ 3.47. 

Another quantity of interest is the average wave number, or spectral average [581 

nana 



iV(t) = v / 1 (9) 

The spectral average is a measure of the number of active modes. In the absence of 
energy cascade, N(t) is limited by the number of linearly unstable modes. 



Figure 12 shows the time-averaged spectral entropy (S) t (top panel) and the time- 
averaged spectral average {N) t (bottom panel) for the kinetic (black circles) and mag- 
netic (red triangles) energy spectra as a function of rj. It can be concluded from the 
top panel that the spatial structures of the velocity field become less complex with the 
onset of the on-off intermittent dynamo at rj ~ 0.053, with a decrease in {Sk) t - This 
is due to the beginning of the action of the Lorenz force on the momentum equation. 
From rj — 0.055 to rj = 0.04 there is a strong decay not only in {Sk) t , but also in (S m ) t , 
reflecting the increase in the frequency of occurrence of large-scale coherent mean- 
field structures in the coherent-incoherent intermittent dynamo as 77 is decreased. The 
spatial complexity of both the velocity and magnetic fields, then, starts to increase 
from rj ~ 0.04 to rj ~ 0.01. This tendency coincides with an increase in the magnetic 
spectral average {N m ) t (bottom panel), which reflects a rise in magnetic energy in 
smaller scales. Note that the kinetic spectral average {Nk) t is always close to the 
energy injection scale k = 5, and the magnetic spectral average is lower due to the 



high energy peak in k = 1 (see Figs. 10 and 11 1 
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Figure 11: Peak height of modes k\ (dashed lines) and fcs (solid lines) of the 
time-averaged magnetic energy spectra as a function of i] in (a) linear and (b) 
linear-log scales. 
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Figure 12: (Colour online) (a) Time-averaged spectral entropy as a function of 
77 for the kinetic (black circles) and magnetic (red triangles) energy spectra; (b) 
time-averaged spectral average as a function of 77 for the kinetic (black circles) 
and magnetic (red triangles) energy spectra. 
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The spectral entropy S m (t) is able to distinguish the two types of intermittency 
displayed in Fig. [5] of section 3.3 Figure 13 a) shows the time series of the spectral 
entropy S m (t) and magnetic energy (-B 2 ) /2 for the on-off intermittency at rj = 0.053. 
There is little correlation between the two S6ri6s, sinc6 Sm 

(t) oscillates erratically even 

during phases of grand minima in the magnetic energy time series. Figure [ljjjb) plots 



the same time series as Fig. 13 a), but for the coherent-incoherent intermittency at 
r\ = 0.05. There is now strong correlation between S m (t) and (-B 2 ) /2, as the entropy 
decreases whenever the magnetic energy increases abruptly during spatially coherent 
bursts, indicating a high-variability not only in the energy amplitude, but also in the 
spatial complexity. 



4 Conclusions 

We have described an intermittent route to mean-field dynamos as a function of the 
(low) magnetic Prandtl number Pr. The main contributions of this work are the 
quantification of spatial complexity of the magnetic field patterns from 3-D compress- 



ible MHD dynamo simulations (Figs. 12 and 13 1, and the identification of a new 



type of intermittency characterized by the switching between coherent and incoherent 
large-scale structures [Figs. [5^b),[7] and 13 b)]. The mechanisms responsible for this 



intermittency have not yet been explored and, in future works, we expect that a bet- 
ter description can be achieved by means of the analysis of phase synchronization of 
Fourier modes [28] [15] • 

This paper is also an addition to a series of previous works on the onset of nonlinear 
dynamo action in ABC flows. Kinematic and nonlinear dynamos were studied by 
Galanti et al. [18] for several choices of the forcing scale kf, and Pr ranging from 1 to 
13, with low kinetic Reynolds number Re (up to 20); A bifurcation study with low-J?e, 
Pr — 1, and kf = 1 was conducted by Seehafer et al. |50| ; Sweet et al. [551156] detected 
a blowout bifurcation responsible for on-off intermittency for kf — 1, Re = 6.3 and 
Pr close to one; Brummell et al. [9] studied kinematic and nonlinear regimes with Re 
between 50 and 100 and kf = 1 in a time-dependent ABC flow previously introduced 
by Galloway & Proctor [20]; Mininni [32] investigated the inverse energy cascade at 
small Pr (down to 0.005) and kf — 3, with Re varying from 11 to 6200; Alexakis & 
Ponty Q] studied the effect of the Lorenz force on on-off intermittency in ABC flows 
with kf = 1 and varying both Re and Rm. In contrast to all the aforementioned 
papers which performed incompressible MHD simulations, in this paper we opted 
for a compressible code. It is not yet clear whether the compressibility is a crucial 
feature of the dynamics observed, though certainly compressibility is a feature shared 
by the solar dynamo! Our simulations were performed with kf — 5, Re ~ 100 (or 
Re ~ 500 for the box-scale Re) and Pr g [0.0625,0.5]. Despite the use of a model 
with simple geometry, reasonably small kinetic and magnetic Reynolds numbers, and 
without differential rotation, our results reveal some qualitative resemblance to certain 
aspects of the solar dynamo and we expect them to be useful in the analysis of more 
realistic models. 
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